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Abstract. We analyze bimolecular reactions that proceed by a long-ranged reactive interaction, using a 
field theoretic approach that takes into account fluctuations. We consider both the one-species, A + A — > 
reaction and the two-species, A + B — » reaction. We consider both mobile and immobile reactants, both 
in the presence and in the absence of adsorption. 

PACS. 82.20.Mj Nonequilibrium kinetics - 05.40.+J Fluctuation phenomena, random processes, and Brow- 
nian motion - 82.20.Db Statistical theories (including transition state) 



1 Introduction 

Most chemical reactions occur through a local mechanism, 
where reaction occurs at a finite rate only when reac- 
tants are closer than a capture radius. In some unusual 
cases, however, the reaction can proceed by a long-range 
mechanism Q . Exciton decay by a multipolar interaction 
is one such example. In these long-range cases, the rate 
of reaction between two reactants is distance dependent: 
w(r) = 7/r" for two reactants separated by a distance r. 
Here the strength of the interaction is given by 7 = Rq / r, 
where Rq is the Forster radius, and r is the decay life- 
time. A variety of electronic energy transfer reactions can 
be modeled by this form [||. More generally, an exponen- 
tial form of the reaction rate may be considered, as in 
electron-hole recombination due to wave function overlap 
or barrier tunneling. The exact form of the interaction 
can be quite complex and so even these shorter-range 
cases may equally well be modeled by the multipolar form, 
e.g. with n — 12 j|J|. A wide variety of physical systems 
exhibit such long-range reaction, including biological elec- 
tron transport, fluorescent decay of electron-hole pairs in 
amorphous semiconductors, scavenging of trapped elec- 
trons in organic glasses, and decay of localized electronic 
states is mixed organic solids. 

The case of long-range reaction has received less theo- 
retical attention than has the more common case of local 
reaction. The one-species, A + A — > reaction between 
immobile reactants without adsorption has been investi- 
gated by placing bounds on a mean-field treatment B. 
An interesting power-law decay of the concentration was 
found in d spatial dimensions at long times: c(t) ~ at~ d l n 
as t — > 00 for n > d. An extension to the A + B — > 
case gave the scaling c(t) ~ a t~ d /( 2n ~ d ) as t — > 00 for 
n > d. The single-species reaction, A + A — > reaction 
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between immobile reactants in the presence of adsorption 
has been investigated with a similar treatment |Q . Adsorp- 
tion corresponds to the creation of reactants, — > A, with 
the conventional adsorption rate J. A non-trivial scaling 
of the concentration with the adsorption rate was found: 
c(J) ~ a J d /i n + d ) as J — > 0, n > d. In principle, however, 
these mean-field results could be modified by renormal- 
ization effects in low spatial dimensions. 

In this paper, we apply the renormalization group ap- 
proach to reaction kinetics to derive the asymptotically 
exact behavior of these long range reactions. There are 
eight distinct physical cases, considering all possibilities 
of one- or two-species reactions, with or without adsorp- 
tion, and mobile or immobile reactants. Our paper is orga- 
nized as follows. In section we review the field theoretic 
approach, displaying the actions appropriate for both the 
one- and two-species reactions in the general case. In sec- 
tion ||, we derive the behavior of mobile reactants in the 
long-time limit in the absence of adsorption. In section 
we derive the long-time behavior of immobile reactants 
in the absence of adsorption. In section we derive the 
steady-state behavior of mobile reactants in the limit of 
small adsorption rates. Finally, in section ^ we derive the 
steady-state behavior of immobile reactants in the limit 
of small adsorption rates. We summarize our results in 
section |?]. 



2 Field Theoretic Formulation of Reaction 
Dynamics 

We consider the general case of a bimolecular reaction in 
d spatial dimensions. The reaction occurs as a result of a 
long-range interaction. For convenience, we define the re- 
action to proceed on a cubic lattice. We use a continuous- 
time master equation to define how the probability of any 
given configuration of reactants changes with time. The 
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master equation for the A + A — > reaction is 

dP({m},t) _ 



dt 



D 



"£[(n j + l)P(. 



l,n j + l,...,t)-n i P} 



^y^ y w ik [{ni + l)(n k + 1)P(. ..,n,i + l,n k + l,...,t) 



- n 4 n fe P] 



(i) 



where h is the lattice spacing, J is the rate of adsorp- 
tion, D is the diffusivity, rii is the number of A species 
on lattice site i, and Wi k the distance-dependent reactive 
interaction. The summations over i and k are over all sites 
on the lattice, and the summation over j is over all near- 
est neighbors of site i. For simplicity, we choose to place 
the species initially at random on the lattice, with average 
number density uq. 

Using the coherent state representation, we map this 
master equation onto a field theory We find that 

the reactant concentration, averaged over the random ini- 
tial conditions and the random statistics of the reaction- 
diffusion process, is 



CA(x,t) = (a(x,t)) , 



(2) 



where the average is taken with respect to exp(— Saa)- 
The action is given by 

Saa = J dd * J o dt s ( x ' *) i dt ~ Dv2 + 5 W] a(x ' f ) 
-no J d d xo(x,0)- J J d d x J dta(x,t) 

+i J d d xd d y J dta(x,t)o(x,*)ii;(|x-y|)a(y,t)a(y,t) 

+ J d d ^d d y J dta(-x,t)a(-x,t)w(\x.-y\)a(y,t) . (3) 

The time U must be larger than the longest time for which 
we wish to make predictions. Special cases of this general 
formulation are given as limiting values. For immobile re- 
actants, for example, we set D = 0. If there is no adsorp- 
tion, we set J = 0. 

For dissimilar reactants, the formulation is slightly dif- 
ferent. The master equation is given by 

aP({mO,{ni},i) 



dt 



D 



y^\( m i + 1 )- p (- • • > m i - l,rnj + 1, .. .,*) - rriiP] 



D 



+ '^2wik[(m i + l)(n k + 1)P(. . . ,m,i + l,n k + 1, . . . ,t) 



- rriin k P] 
-Jh d Y J [P{---,m i -l,...)-P] 

i 

-jh d j2[ p (---^i-i,...)-p] , 



(4) 



where now rrii and rii are the number of A and B species, 
respectively, on lattice site i. We have assumed for sim- 
plicity that the diffusion constants of the two species are 
the same. So as to reach an interesting scaling regime, we 
have set the adsorption rates of the two species to be the 
same. Finally, we will set the average initial concentrations 
of the two reactants to be the same. 

We derive the field theory from this master equation 
using the coherent state representation ]Tl|,[l2|. The aver- 
aged concentrations are given by averages over two dis- 
tinct fields: 



CA(x,i) = (o(x,t)) 
cs(x,t) = (&(x,t)) , 



(5) 



where now the average is with respect to exp(— Sab)- The 
action for dissimilar species is given by 



5' 



AB 



d d yi I dt a(x, t) [d t - DV 2 + 8{t)] a(x, t) 



d d x / dt 6(x, t) [d t - DV 2 + 5(t)] 6(x, t) 
Jo 

n J d d x [a(x, 0) + 6(x, 0] 

J J d d x J dt [a(x, t) + 6(x, t)] 



d d xd d y / dto(x,t)a(x,t)ti;(|x- y |)6(y,t)6(y,i) 
Jo 

J d d xd d y j dto(x,t)a(x,t)«;(|x- y |)6(y,t) 



+ Jd d xd d yJ^ rft6(x,t)6(x,fH|x-y|)a(y,t) . (6) 

The long-ranged interaction shows up in the non-quadratic 
terms of this field theory. The vertex that we must con- 
sider for the single-species reaction is shown in figure [l]. 
Similarly, the vertex we must consider in the two-species 
reaction is shown in figure |. We will use both a regular- 
ized interaction 



w R {r)= 1 /{r 2 + R 2 Y' 2 



(7) 



and a pure power-law form of the interaction. The Fourier 
transform of the interaction is given by 



w)^(k) = J d d r exp(ik • r)wR(r 

_ 7 7r d / 2 (fci?/2)("- d )/ 2 



ik 



R n - d r(n/2) 



2K {n _ d)/2 (kR) , (8) 



Jeong-Man Park and Michael W. Deem: Field Theoretic Approach to Long Range Reactions 



3 




Fig. 1. The vertex in the field theory of the A + A — > 
reaction. Time increases in the direction of the arrows. The 
wavy line indicates reaction with the distance-dependent rate 
w(r). 




Fig. 2. 

reaction. 



The vertex in the field theory of the A + B 



where F(x) is the Gamma function, and Ki(x) is the mod- 
ified Bessel function of the second kind. At long-times, for 
low adsorption rates, the reactants tend to be widely sep- 
arated, and the regularization does not affect the reaction 
kinetics. In the absence of regularization, we find 



lim [w R (0) - tSfl(fc)] - w(0) 

R — *0 



w(k) 

~fir d / 2 r[{d- 



jn—d 



n)/2] k- d . (9) 



r(n/2) 

This equation is well-defined by the Fourier integral for 
d < n < d + 2. It is derived for all other n by analytic 
continuation from the relation FT{V 2 w} = —k 2 FT{w}. 



3 Mobile Reactants without Adsorption 

In this section we address the simplest case of bimolecu- 
lar reactions between mobile reactants in the absence of 
adsorption. That is, we set J = in the field-theoretic 
actions. We use renormalization group theory to calculate 
the asymptotic reactant concentrations in the long-time 
limit. 





a) 



b) 



Fig. 3. The local vertex and vertex correction in the field 
theory of the A + A — » reaction, a) The effective local inter- 
action, into which the non-local vertex flows, b) The vanishing 
correction to the non-local interaction. 



3.1 The A + A 



Reaction 



For the single species reaction, we find the flow equations 
to be 



din no 

dl 
dlnD 

dl 
din 7 

dl 



d 

z-2 
z — n 



(10) 



where z is the dynamical exponent. These flow equations 
are exact to all orders. Time renormalizes as t(l) = e t. 
To reach the fixed point, which requires the propagator 
to reach a fixed point as well, we set z = 2. We see that 
the strength of the long-range reaction, 7, appears to be 
irrelevant when n > 2. What this result actually implies 
is that the long-range nature of the reaction is irrelevant. 
Certainly, though, the reaction event itself is relevant. We 
can, therefore, use the effective interaction w(r) — XS(r). 
This replacement is possible only when the interaction 
is integrable, n > d, which we always require. For non- 
integrable interactions, the reactant concentration decays 
to zero immediately. The scaling for the effective reaction 
rate is 

*!E± = z-d--t-. (if) 

dl 4ttD y ' 

This flow equation is a one-loop expansion for small A. 
The flow equations |l0| and [ll] are derived by considering 
the vertex corrections shown in figure |^. 

The renormalization group transformation relates the 
original system at long times and low concentrations to 
another, renormalized system at short times and high con- 
centrations. In fact, we integrate the flow equations until 
the renormalized time is short enough so that we can use 
simple mean field theory JToj] . The matching time 



te 



./'* Z (l)dl 



to 



(12) 



is chosen to be on the order of h 2 /(2D). Calculating the 
renormalized reactant concentration at short times with 
mean field theory, we find 



c[t(i*y,n = 



1 



l/n (l*) + X(l*)t(l*) 
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a* 



b* 




a) 



b) 



Fig. 4. The vertex corrections in the field theory of the A + 
B — » reaction, a) The effective local interaction, into which 
the non-local vertex flows, b) The vanishing correction to the 
non-local interaction. 



1 



X(l*)t(l*) 



as I* 



(13) 



The concentration of the original system is related to that 
of the effective, renormalized system by scaling: 



c(i) 



-dl 



[t(l);l] 



(14) 



Combining equations [l3] and |l4|, and using equation [l2] to 
express the result in terms of t rather than I* , we find 



As in section the prediction in equation [Hj is limited 
to the diffusion-limited regime, which section shows 
occurs for d > 4 or d < 4 and n > 4. 

In summary, for mobile reactants, the long-range na- 
ture of the reaction is irrelevant in the diffusion-limited 
regime. The long-time scaling of the concentration is the 
same as that of a reaction with a local interaction for both 
the one- and two-species reactions. 



4 Immobile Reactants without Adsorption 

The case of immobile reactants in the absence of adsorp- 
tion, D = and J = 0, is rather different. This is be- 
cause for immobile reactants, the long-range and short- 
range cases are quite distinct. This is simply because a 
short-range reaction stops as soon as no more than one 
reactant occupies each lattice site. The long-range reac- 
tion, on the other hand, continues until either zero or one 
reactant remains in the entire system. As we will see, this 
essential difference leads to a more complicated matching 
procedure. 



4.1 The A + A 



Reaction 



At 



c A (t) 



ln(t/tp) 
8nDt 



,d>2 
,d = 2 



as t 



oo. 



(15) 



s -d/2 



Here, we have used the fact that X(l) goes to a fixed point 
value for d < 2. To first order in 2 — d, we find from 
equation |ll] that A* = 4tt£>(2 - d). 

The predicted decay in equation ^| comes from as- 
suming that the reaction will be diffusion limited at l ong 
times. This is generally true, but the results of section 11 
show that when d < 2, the long-range nature of the reac- 
tion is relevant when n < 2. Our prediction [l5|, therefore, 
is limited to the cases of d > 2 or d < 2 and n > 2. 



3.2 The A + B 



Reaction 



For the two species reaction, we find the same flow equa- 
tions as in equations 10 and O. The vertex corrections 
that lead to this result are shown in figure |^. As with the 
single species reaction, the long-range nature is irrelevant 
when n > 2. At long times, the behavior of the long-range 
case is identical to that of the local case. This implies, as 
before, that z = 2. The matching is somewhat more com- 
plicated than that of the single species reaction, due to 
the segregation between the A and B species that occurs. 
The result is HI 



c A (t) = c B {t) 



j_ 

\i 



,d > 4 



V"o 1 <• A ^ A 

^{8TrDt) d / i ' 1 - " ^ ^ 



as t 



oo. 



(16) 



As before, we use renormalization group theory to analyze 
the field theory in the long-time regime. Now, however, 
we must be careful with the definition of the interaction. 
We use the regularized version in equation [?]. Since we 
are interested in the long-time regime, where the particle 
density is low and the particles are widely separated, we 
expand the interaction for large r: 



To 



72 



(17) 



where 70 = 7. We find the flow equations to be 

dh\TlQ 

dlwyi 

= z — n — 1 . 

dl 



(18) 



We immediately see that the higher order terms in the 
reaction rate, 7, for i > 0, are less relevant than 70. So as 
to reach a fixed point, we set z = n. The flow equations 
for no and 70 are then asymptotically exact. 

While the regularization of the interaction w(r) is ir- 
relevant for long-times, it is important in the matching 
limit. Indeed, we integrate the flow equations until the 
density is on the order of c[t(l*);l*} w 1/R d . In that limit, 
wc find 



c(W*) 



1 



1 



no 1 (I* 



w(0)t Q w{0)t o 



as I* 



(19) 



where w(0) = (const)7i? d " is not a function of I. We 
then find 



c(t) 



Wu! (to) 



-d/n 



,n > d 
,n < d 



as t — > 00. 



(20) 
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As indicated, the specific value of the cutoff affects the 
prefactor, but not the exponent, in the long-time scaling. 
When the interaction is not integrable, n < d, the concen- 
tration immediately decays to zero in the limit of a large 
system size. 

While our treatment of this case has required use of a 
cutoff, the reaction process defined by the master equation 
is well-defined even in the absence of a cutoff. This is be- 
cause any reactants experiencing the infinite reaction rate 
w(0) simply react immediately at the onset. Indeed, in the 
limit R — > and h — > 0, there is an exact scaling relation 
that relates the system at long times to a renormalized 
system at short times. If we rescale time as t* — ^n^ d t 



and space as x* 



lid 



x, we find that 
c(x, t) = noc*(x*, t*) 



(21) 



where c* is the exact solution of the master equation [l] for 
no = 1 and 7 = 1. This scaling relation explicitly demon- 
strates that a cutoff is not required. Indeed, this scaling 
relation is exactly what the renormalization group treat- 
ment is intended to derive. Essentially, this problem has so 
few time and length scales that the renormalization group 
approach does not provide any additional simplification. 

On a deeper level, the scaling relation |l| implies that 
there is no non-trivial renormalization of the effective in- 
teraction for this problem. Mean field theory, therefore, 
will be accurate for both long times and short times. If we 
assume that there are no correlations among the reactants, 
we can write 



^ = -c f°° drSy-^cM 
dt J R r" w 



where S d = 2n d / 2 / T(d/2) is the surface area in d di- 
mensions. Integrating this equation using the approxima- 
tion that there is only one reactant per correlated region, 
e(t) w d/[S d R d (t)], we find 



c(t) 



d_ 
S~d 



jnt 
n — d 



-d/r, 



as t 



00, n> d. 



(23) 



This relation gives the same power-law decay as does equa- 
tion |2(| Moreover, it satisfies the exact scaling relation 
This approximate approach has essentially assumed that 
the radial correlation function rises from zero to unity at 
the average distance of interparticle separation. A more 
careful treatment of the correlation function would likely 
lead to a more accurate prefactor for the asymptotic con- 
centration decay, without changing the exponent [Us. 



4.2 The A + B 



Reaction 



For the two-species reaction, we find the same flow equa- 
tions as for the single-species reaction, equation [l8| The 
matching is different, however, because the A and B reac- 
tants segregate in the long-time limit, just as they do in 
the mobile case. The segregation of the reactants requires 
us to adopt a different dynamical exponent. We argue that 



the degree of reaction between a typical reactant and a 
patch of reactants of volume AV separated by a distance 
r is invariant at the fixed point. The typical number of 
reactants in a volume AV is proportional to (AV) 1 / 2 in a 
segregated system, and so we require ^(AV) 1 ^ 2 At/r n to 
be invariant under the spatial and temporal rescaling in 
the renormalization group procedure. This condition gives 
z = n — d/2. 

To analyze the matching regime, we use the mean field 
equations 

dc f 

-^=-c A (x,t) / d d yc B (y,t)w(\x-y\) 
dc B 



dt 



= -c s (x, t) / d d y c A (y, t)w(\x - y|) . (24) 



Since we expect segregation, we rewrite these equations in 
the variables <p — ( C A — c s)/2 and p — (c A + cb)/2. We 
make the approximation of uniform total density: p(x, t) w 
p(t) . This approximation is valid in the limit of weak seg- 
regation, <fi <C p, and in the limit of strong segregation, 
p = \<fi\, except near domain boundaries. With this as- 
sumption, we find 



and 



— = p{<f> * w) - (j>pw(0) 



dp 

— = </>(</) *w) - p 2 w(0) 



(25) 



(26) 



Equation can be solved, since p depends on time only: 



(22) 0(k,i) =0(k J Q)eacp^-[t&(O) -*&(*)] / dt' p{t') 



(27) 

Taking the average of equation|26j and using (4>(r, 0)^(r', 0)) 
n S(r — r')/2, we find 



d(p) _ n 



dt 



)(k) ^cxp j-2[u>(0) - w{k)} J dt'p(t')\^ 



( P 2 )w(0) 



(28) 



where the notation J k stands for J d d k/(27r) d . Note that 
all the terms in equation ^8] scale as e dl I 2 in the matching 
limit. We also know that the initial reactant concentration 
decreases from noe dl to y/n e dl I 2 in a time of the order 
of l/[w(0We dr ] |||. Knowing this, and using the scaling 
relation [fi], we see that the concentrations must scale as 
ca — cb ~ (const)i~ d /( 2n ~ d ) in the long-time limit. 

The exact scaling relation |2l] also applies in this case, 
however, and so there are no non-mean-field effects intro- 
duced by renormalization. Knowingthis, we can analyze 
the approximate mean-field theory [28| more thoroughly to 
obtain an estimate of the prefactor in the long-time de- 
cay. We perform the calculation with the non-regularized 
interaction of equation ^. First, we calculate the average 
fluctuations of <b: 



n Q S d 



nd/o 



2< 



2a f* dt'p{t') 



d/C 



(29) 
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where we have defined £ 



d and 



5.1 The A + A 



Reaction 



2«r(n/2) 



(30) 



Balancing the terms in equation |26|, we find that the two 
right-hand terms must be equal, p(t) ~ |<H0l- Indeed, 
setting <p — at~ 5 , we find § = d/ (2n — d). We also find an 
explicit form for a from equation by noting that since 
is a Gaussian field, (\<f>(t)\) = \2[4> 2 ) /ir] 1 ' 2 . We finally 
obtain 



(P) 



n {n-d)/{2n-d) o -d/{2n-d) 



s d r(d/0] in - d)/i2n - d) 



2< 



x 



d/(2n-a 



. 2n - d , 
as t — > oo, n > d 



1/2 



-d/(2n-d) 



(31) 



This exponent is in agreement with the scaling argument 
of Burlatsky and Chernoutsas |6) . This result can be non- 
dimensionalized exactly as predicted by equation |2lj. Our 
approach here has made use of a self-consistent approx- 
imation. In particular, we integrated equation 25 by as- 
suming that pit) was a given function. We then, however, 
made the identification p = \<j>\ to obtain the prefactor in 
equation [u]. 

As in the single-species case, we can also take a correla- 
tion function approach to the calculation of the long-time 
concentration profile. We have 



dcA 
dt 



= -CA 



dr-Sd^^csM . (32) 



Again assuming that the radial distribution function rises 
from zero to unity at the average interparticle separation 
r = R, we find c A (t) = {n d/[S d R d (t)]} 1 l 2 in the segre- 
gated limit. Integrating equation B3, we find 



CA(t) = c B (t) ~ n Q 



2n 



n — d 

-d/(2n-d) 



-d/(2n-d) 



as t 



n/(2n-d) 



s dj 

n > d. 



(33) 



This result satisfies the exact scaling relation |2l]. As in the 
single-species case, a more accurate correlation function 
approach would likely lead to a better prefactor, without 
modification of the exponent jl4) . 



5 Mobile Reactants with Adsorption 

We now turn to consider the case of long-range reaction of 
mobile reactants in the presence of adsorption. That is, we 
consider the general case defined by the master equations 
|l| or ||, with fl / and J ^ 0. To access an interesting 
scaling regime, we consider the limit J — > and seek to 
understand how the concentration scales with J in this 
limit. 



We find the flow equations for the single-species reaction 
to be 

din no 



dl 
din J 

dl 
dlnD 

dl 
dhij 

dl 



= d 



(34) 



These flow equations are exact to all orders. We see that, 
as for the case of no adsorption, the long-range nature 
of the reaction is irrelevant when n > 2. Of course, the 
presence of the reaction is important, so we set the inter- 
action to the effective value w(r) = A<5(r). We require the 
interaction to be integrable, n > d, as always. The flow 
equation for A is the same as in equation [n]. To reach the 
fixed point, we set z — 2. We integrate these flow equa- 
tions until J(P) = e ( 2 + <i ) r J = J ~ D/h 2+d . For this 
renormalized value of J, the average density is finite, and 
we can use mean field theory. The mean field theory is 
simple in this case, predicting 



c[J(l*),t=oc;l*] = [J(r)/A(P)] 1/2 



(35) 



We match the mean field result to the observed value using 
the scaling relation 

c{J) = e- dl 'c[J{r)-r] . (36) 

We find for the concentration in the limit J — ► 0: 

,d> 2 



:f)" 2 



■/In (Jo/. 7) 
16nD 



1/2 



2 as J -> 0. (37) 



c(J) 



Here, we have used the fact that A(/) goes to a fixed point 
value for d < 2, just as it does in the case of no adsorption. 
These predictions are consistent with those of Rey and 
Droz pq |. Our argument that the long-range nature of 
the reaction is irrelevant is, again, v alid when d > 2 or 
d < 2 and n > 2, as shown in section 3.1. 



5.2 The A + B 



Reaction 



For the two-species reaction, we find the same flow equa- 
tions as for the single-species reaction. We also find that 
the local interaction dominates over the long-range com- 
ponent. The matching, however, is somewhat different. 
Again, this difference is due to the reactant segregation 
that can occur in low dimensions. The appropriate mean- 
field equations in the matching limit are 

Oca 
dt 

^ = DV 2 c B - \c A c B + J B , (38) 



DV c A - Xc A c B + J a 
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where J a and Jb are random, Poisson adsorption fluxes of 
A and B reactants, respectively. Note that (J a) = (Jb) = 
J. The initial conditions are irrelevant in the long time 
limit, since the steady state behavior is controlled by the 
adsorption rate. We assume for simplicity that n = 0. 
These equations can be simplified and analyzed in a fash- 
ion similar to the case of no ads orption [ pr||l^ | . Using the 
same variables as in section 4.2, we find 



dt 



D\7 A 



SJ 



and 



dp 2 2 /2i , J A + Jb 



(39) 



(40) 



where 5 J — {J a — Jb)/^- We can solve equation ^9|to find 

<Kk,t) = / dt'e' D ^ k2 SJ(k,t') . 
Jo 

Averaging equation we find 

= -X[(p 2 )-(^)}+J . 
From equation Ed], we find 



(41) 



(42) 



lim 

t — >-oo 



= lim — 

t—>oo 2 

_ ±f 

2D Jd , 



dt' 



D -4dDt/h 2 
h 2 



(43) 



where fd is a dimensionless numerical constant depending 
on the dimension. Here Io(x) is the modified Bessel func- 
tion of the first kind. For two or fewer dimensions, fd = oo. 
We can bound the solution of equation ^ First, 
we note that since the concentrations are non-negative, 
p > \4>\, and so (p 2 ) > (</> 2 ). Using this, we find 



\) 2 <((p-\4>\) 2 ) 

= (p 2 ) + {4> 2 )-2{p\4>\) 

< (p 2 ) (4> 2 ) 
= j/x . 

Combining the lower and upper bounds, wc find 



(44) 



Jfd 
irD 



1/2 



<(P)< 



Jfd J 

2D A 



1/2 



(45) 



We can immediately conclude that for a local reaction 



(P) = 



d<2 



(46) 



This occurs because any given lattice site is dominated 
either by A or B species due to segregation. The random 
rate of adsorption of the dominant species is so large for 
d < 2 that it swamps the relaxation by diffusion of the 
species to nearest neighbor lattice sites, which is the only 
mechanism for reaction in a segregated system. Since the 



local reaction is ineffective for d < 2, we must reconsider 
our argument that the long-range component of the inter- 
action is irrelevant. Indeed, the results of section |6 . 2j imply 
that equation ^ holds only for n > 2d. 

In the matching limit, for d > z, the flow equations 
lead to 

D(l) = De^ 2 ^ 1 
J(l) = j e {z+d » , 



(47) 



where we set z = 2 to obtain a fixed point for the diffu- 
sivity. We define e = J(l)fd/[2D(l)] and note that e <C 
J(l)/X(l). We see that 



([c A (l)-c B (l)} 2 



4e 



and 



([c a (1)+cb(1)} 2 ) = 



4J(/) 
A(0 



4r 



(48) 



(49) 



We note that equation 
field, 



implies, since is a Gaussian 
.-|V2 



(\c A (l) - c B (l)\) 



(50) 



This implies that ca(1) and cg(i) fluctuate in the same 
way: c A (l) ~ c B (l) + 0(y/e). Since (p 2 (l)) = (0 2 (O> + 
J(/)/A(0, we find (c A (l) 2 ) = (c B (l) 2 ) = J(l)/X(l) + O(e). 
Using (a(l) 2 } = (ci(l)} 2 + ([Sci(l)] 2 } and equations |4^ and 
H we find ([Sc A {l)] 2 ) = ([Sc B (l)} 2 ) = 0(e). Since the fluc- 
tuations are much smaller than the mean, we can finally 
conclude that 



(ca(1)} = (c B (l)) 



J(l) 



1/2 



as 0, d> 2 . (51) 



Since the scaling predicted by the flow equations is trivial, 
this result holds for the observable concentrations as well: 



(ca) 




,d>2 
,d<2 



as J — > 0. 



(52) 



6 Immobile Reactants with Adsorption 

We now consider the final case of immobile reactants in 
the presence of adsorption. We consider the regime J — * 0, 
which leads to non-trivial scaling of the reactant concen- 
tration with the adsorption rate. As in section |^, we will 
find that the decay rate depends sensitively on the inter- 
action w(r). 



6.1 The A + A 



Reaction 



For the single species reaction, we find the sa me fl ow equa- 
tions and dynamical exponent as in section 4.1, with the 
additional equation for the adsorption rate 



d In J 
dl 



z + d 



(53) 
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As in section 4.1, we find that while simple scaling ar- 
guments predict the local reaction mechanism to domi- 
nate over the long-range reaction mechanism, this is actu- 
ally not the case. In the low-concentration, small- J limit, 
the reaction proceeds only by the long-range mechanism 
for immobile reactants. In the matching limit, we find 
c[J(l*);l*] = [J(P)/J ] 1/2 , where J w jh~ d - n . The ob- 
servable concentration is given by the scaling relation c( J) = 
e~ dl * c[J(l*); I*}, where I* is chosen so that J (I*) — Jo- 
Performing the matching, we find 



In the long-time limit, the average density goes to a con- 
stant. In this limit, then, 



c(J) 




N d/(n+d) 

J oJ as J 

3 ,n < d 



0. 



(54) 



This result is the same as that found by Oshanin and 
coworkers using an approximate, mean-field approach @. 

We can also find an exact scaling relation for this prob- 
lem. Defining x* = ( J/ 7 ) 1 /( Tl+d )x and t* = J{j/ J) d ^ n+d h, 
we find 

/ t\ d/(n+d) 

c(x,t)=(-J c*(x*,r), (55) 

where c* satisfies the exact master equation |l| with J = 
1 and 7 = 1. This scaling relation implies that (c) = 
(const) (j/ 7 )<V(«+d) = (const)(J/J ) d/(n+d) , which is ex- 
actly what we found in equation B4l 



6.2 The A + B 



Reaction 



For the two-species reaction, we find the sa me fl ow equa- 
tions and dynamical exponent as in section 4.2, with the 



addition of equation |53|. In the matching limit, we can use 
a similar type of mean field theory: 



— = p(4> * w) — <fipw(0) + 6 J 
at 



(56) 



and 



dp 
dt 



*w)-p 2 w(0) + (J A + J B )/2. (57) 



Equation pq can be solved, since p depends on time only: 



dt"[w(0)-w(k)}p(t") 



4>(k, t) = / dt' exp 
Jo 

xSJ(k,t') 



(58) 



Taking the average of equation |57j and using equation [38 
and (SJ(k,t')SJ{k',t")) = (J/2)5(t - t")(2ir) d 5(k - k'), 
we find 



M = Jf t d , 
dt 



xw{k) (exp (-2^ dt"[w(0) - w(k)]p(t")\ 
- (p 2 ) w(Q) + J . (59) 



= 



J 



)(k) 



— p 2 w(0) + J 



(60) 



, k 2[w(0) - w(k)]p 

Using the explicit expression ^ for the interaction, we can 
evaluate the first integral: 



w(k) 



A 



= (const) / dkw{k)k d -^ , 
k w{0) - w(k) J 

< oo, d>C- (61) 

Here A = 2n/h is the cutoff in Fourier space, and we have 
concentrated on the small-fc, large-distance properties of 
this integral. Since £ = n — d, we see that the integral con- 
verges to a finite value only for n < 2d. When the integral 
converges, the first term in equation |6^ is subdominant, 
because p is large in the matching limit, and by perform- 
ing the matching we find the expected mean-field behav- 
ior. In the case of n > 2d, the average density diverges, 
as we can see explicitly from the bound (p) > (\(f>\), since 



result is 



ca = c B 



(const) L dkk d ^ 1 diverges for n > 2d. The final 



3C 



d/(n+d) 







,n > 2d 

,d<n<2d as 0. (62) 

, n < d 



This result satisfies the exact scaling relation |55|, as it 
must. 

Interestingly, if the interaction decays too quickly with 
distance, n > 2d, then segregation occurs, and the adsorp- 
tion swamps the reaction rate. This effect can only occur 
in the A + B — > case, where the density of a majority 
species can grow without bound on a lattice site. If the 
interaction is sufficiently long-ranged, d < n < 2d, then 
segregation is prevented, and the reaction proceeds as in 
the homogeneous, single-species case. As always, if the in- 
teraction is not integrable, n < d, the concentration is 
forced to vanish in the limit of a large system. 



7 Summary 

Using a field-theoretic formulation of reaction kinetics, we 
have derived the asymptotic behavior of long-ranged, bi- 
molecular reactions. When the reactants are mobile, the 
long-ranged nature of the reaction is usually irrelevant. 
Conversely, when the reactants are immobile, the long- 
ranged reactive interaction is crucial for the reaction to 
proceed. In the A + A — > reaction, concentration fluc- 
tuations are only mildly important. The reactants in the 
A + B — > reaction, however, typically segregate in the 
absence of adsorption. Adsorption tends either to homog- 
enize or to drive the system to infinite density in the 
A + B — > reaction. 
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